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Abstract 

In this paper, we further investigate the global dynamics of a stochastic 
differential equation SIS (Susceptible-Infected-Susceptible) epidemic model re¬ 
cently proposed in [A. Gray et ah, SIAM. J. Appl. Math., 71 (2011), 876-902]. 
We present a stochastic threshold theorem in term of a stochastic basic repro¬ 
duction number Rq : the disease dies out with probability one if Rq < 1, and 
the disease is recurrent if Rq ^ 1. We prove the existence and global asymptotic 
stability of a unique invariant density for the Fokker-Planck equation associ¬ 
ated with the SDE SIS model when Rq > 1. In term of the profile of the 
invariant density, we define a persistence basic reproduction number Rq and 
give a persistence threshold theorem: the disease dies out with large probabil¬ 
ity if Rq ^ 1, while persists with large probability if Rq > 1. Comparing the 
stochastic disease prevalence with the deterministic disease prevalence^ we dis¬ 
cover that the stochastic prevalence is bigger than the deterministic prevalence 
if the deterministic basic reproduction number Rq > 2. This shows that noise 
may increase severity of disease. Finally, we study the asymptotic dynamics of 
the stochastic SIS model as the noise vanishes and establish a sharp connection 
with the threshold dynamics of the deterministic SIS model in term of a Limit 
Stochastic Threshold Theorem. 
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1 Introduction 


Transmission of diseases are usually described in compartmental models 6 9,12,19 . 


Most of the early work on mathematical epidemiology focused on deterministic models 
However, stochastic noise plays an indispensable role in transmission of 


19 


diseases, especially in a small total population. Therefore it seems more practical to 
consider stochastic epidemic models. 

One way to describe the compartmental stochastic epidemic models is via stochas¬ 
tic differential equations (SDE), which change parameters involved in deterministic 
models randomly. This is the so-called parameter perturbation method of deriving an 
SDE model from its deterministic counterpart [^[8,18 . Such environmental noise re¬ 
flected in specihc parameters of the stochastic epidemic models may indicate whether 
severity of disease increases or decreases based on what the deterministic compartmen¬ 
tal epidemic models predict and thus influence of these parameters is well captured. 

It is well known that epidemic threshold theorem holds for most deterministic 


compartmental epidemic models 12 19 : the disease free equilibrium (DEE) is globally 
asymptotically stable if the basic reproduction number i?o ^ 1, and an endemic 
equilibrium exists and is globally asymptotically stable if i?o > 1- However, there 
seems no global threshold result for SDE epidemic models in literature. 


In 18 , a stochastic differential equation SIR (susceptible-infected-recovered) model 


is discussed and a sufficient but not necessary condition for local almost sure asymp¬ 


totic stability of DEE is presented. Later, in 16 , an improved sufficient condition is 
given by considering an SIRS model which specializes the SIR model in 18 . Both 


of these results are derived by constructing Lyapunov functions. In j^, a sufficient 
and necessary condition for local almost sure asymptotic stability of DEE is proved 
by considering the linearized system. 


Perturbing the constant fd in the deterministic SIS type model 


( 1 . 1 ) 


dS{t) 

dt 

dm 

dt 


=fiN - (3S{t)I{t) + ■ylit) - fiS{t), 

=(3S{t)I{t)-{fi + mit), 


randomly by /3dt = /3dt + adB{t), Gray et al. proposed the following SDE SIS 
model 


( 1 . 2 ) 


dS{t) =[/riV - (3S{t)I{t) + 7J(f) - pR(f)]df - aS{t)I{t)dB{t), 
dl{t) =[/3S{t)I{t) -{fi + 'j)I{t)]dt + aS{t)I{t)dB{t) 


with initial values So+Iq = N. Over a long period of time, the average total population 
size is assumed to be constant N. Here S(t) and I(t) are the number of susceptibles 
and the number of infectives at time t, respectively; B{t) is a Brownian motion; p 
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is the per capita (birth) death rate, 7 the recovery rate, {3 the disease transmission 
coefficient, and a the standard deviation of the noise. For more details of the model, 
we refer the reader to [^|^. 

Since S{t) + I{t) = N, which is a constant, it reduces to studying the following 
one-dimensional SDE: 

(1.3) d/(f) = I{t) ^ ^ + a[N - I{t)]dB{t)^. 


In the global existence, uniqueness, boundedness and positiveness of (1.3) are 

proved, and the counterpart of the basic reproduction number in SDE SIS model 

BN a 6 

is dehned by := -—-r. When Rq < 1 and cr^ ^ —, or cr^ > 


max < — 




/i-Fy 2(p-F7)’ 

, it is shown that the disease dies out with probability one. When 


iV’2(p + 7) 

Rl > 1 , it is demonstrated that the disease persists in the sense that 


(1.4) 


liminf J(t) ^ /*(cr) ^ limsup/(f), 

t—^00 




where R{a) = ^ ( y + d) ~ {P ~ o'^N)) ■ Moreover, a unique stationary 


distribution of model (1.3) is proved to exist and formulae for the mean and variance 
of the stationary distribution are presented. Numerical simulations reveal that as 
long as Rq < 1, the disease will die out in the long run. Thus a conjecture on the 
extinction of the disease is proposed in 1 ^: 


Conjecture 1 . 1 . If Rq < 1 and 

with probability one. 


jB 

— -r ^ > —, then the disease will die out 

2(/i + 7) iV 


In this paper, we further investigate the dynamics of the stochastic SIS model 


(1.2). Instead of using the standard approach of constructing a Lyapunov function, 
we use Feller’s test for explosions (see Lemma 6.2 in the Appendix) to establish a 
global threshold theorem for model ( |1.2 ). We prove that the disease dies out with 
probability one if Rq < 1 , and the disease is recurrent in the sense that the process 


For a precise dehnition of recurrence 

Although the 


of disease transmission is recurrent if Rq ^ 

in this paper, we refer the reader to Lemma 6 ^ in the Appendix, 
existence of the stationary distribution is presented in [^, the prohle of the station¬ 
ary distribution, which contains more information about the SDE SIS model, is not 
addressed. Instead of studying the SDE, we investigate the Fokker-Planck equa¬ 
tion (FPE) associated with (1.3) and solve the invariant density. By constructing a 
proper Lyapunov function, we show that the invariant density is unique and globally 
asymptotically stable. We further dehne a persistence basic reproduction number Rq, 
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disease-extinction with large probability and disease-persistence with large probability 
in term of the profile of the invariant density. We show that the invariant density 
peaks at 0 if Rq < 1, while peaks at some positive number R{a) if Rq > 1. Com¬ 
paring the stochastic disease prevalence with the deterministic disease prevalence, we 
discover that the stochastic prevalence is greater than the deterministic prevalence 
if the deterministic basic reproduction number Rq > 2. These results reveal that 
stochasticity favors disease extinction if the basic reproduction number is small, but 
enhances severity of disease if the basic reproduction number is large. By analyz¬ 
ing the asymptotic behavior of the invariant density as a —)■ 0, we establish a sharp 


connection with model (1.1) via a Limit Stochastic Threshold Theorem. 


This paper is organized as follows. In the next section, we present the global 
stochastic threshold theorem. We derive the existence, uniqueness, global stability 


and an explicit formula of an invariant density of the FPE associated with (1.3) in 


Section 3. Then we dehne the persistence basic reproduction number and give a per¬ 
sistence threshold theorem in Section 4. Finally, we study the asymptotic behavior of 
the invariant density as the noise vanishes and establish the Limit Stochastic Thresh¬ 
old Theorem in Section 5. We list some preliminary results on Feller’s test and FPE 
in the Appendix. 


2 Stochastic Threshold Theorem 


In this section, we prove a stochastic threshold theorem: the disease dies out with 
probability one if R^ < 1, but is recurrent in the sense that the process of disease 
transmission is recurrent if Rq ^ 1. 

For the reader’s convenience, we hrst restate the global existence of a unique 


bounded positive solution of (1.3) 


Lemma 2.1. For any given initial value /(O) = /q € (0, A^), the SDE (1.3) has a 
unigue global positive solution I{t) G (0, N) for allt ^ 0 with probability one, namely, 

F{I{t) e (0,A^) : V f ^ 0} = 1. 

Now we state the following stochastic threshold theorem, which in particular, 
verihes Conjecture 1.1 in Section 

Theorem 2.2 (Stochastic Threshold Theorem). 

(1) If Rq < 1, for any given initial value /(O) = /q G (0,A^), 

p| lim I{t) = o| = 1. 

L t^oo j 

In other words, the disease dies out with probability one. 
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(2) If Rq ^ 1, for any given initial value /(O) = Iq E (0,iV), 


P<^ sup I(t) = iV ^ = P<^ inf I(t) = 0^ = 1. 
lo^t<oo J loa<oo J 

In particular, the process R is recurrent: for every x G (0,iV), we have 

P{/(f) = X ■. 3 f G [0, cxd)} = 1. 

In other words, the disease is recurrent. 


Proof. By Lemma 2.1, let Y(t) = g{I{t)), where (7 : (0, iV) —)■ M is defined by 

e 


( 2 . 1 ) 


giO = log 


iv-e 

Using Ito’s formula, it is easy to verify that Y(t) solves the following SDE 


( 2 . 2 ) dU(f) = 


(^flN - /i - 7 - -cr^A^^ j - {t^ + 7)e’^l‘l + 


1 + 


df + aNdB{t). 


In fact, by (|1.3|), 
dY{t) = 




dt + a{I{t))g'{I{t))dB{t), 


where a(0 = (t{N - OC, &(0 = - /i - 7 )^ - h'iO = + t and 


Thus 


iv-e e 


(2.3) dU(f) = 


[on -g-T- + X-NIiV 


N - I(t) 


dt + aNdB{t). 


_ _ 

Substituting I{t) = g~^{Y{t)) = -- 77 ^^ into (2.3), we obtain (2.2) 


1 + e^h) 


Next, by Lemma 2U and Lemma 6 N in the Appendix, we have the global existence 
of a unique solution to ( 2 . 2 ): 


Theorem 2.3. For any given initial value U(0) = lo ^ 1^, the SDE (2.2) has a 
unigue global solution Y{t) G M for all t ^ 0 with probability one, namely. 


(2.4) 


P{U(t) G M : V t ^ 0} = 1. 
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By Theorem 2.3, the process Yt is conservative (for its dehnition, see p.l53 in [lO]). 


The scale function dehned before Lemma 6.2 in the Appendix for (2.2) is given 


by 


i){x) = / 


with 

0(0 = exp 

= exp 


2 




(^(3N - - (h + iV + 


1 + e^ 


(/3iV -/i - 7 - e + (e^ - 1) - 2 log 


dr 

e« + 1 


It is obvious that 'ip{oo) = oo. Note that 


0(0 ~ exp I - - 

/3N 


2(/i + 7) 


+ 21 og 2 >, as — cxD. 


Recall that = 

/i + 7 2(/i + 7) 


. If Rq < 1, then fiN — /i — 7 — -cr^A^^ < 0, which 


implies that 'ijj{—oo) > — 00 ; if Rq ^ 1 , then (]N — /x — 7 — ^ 0 , which implies 

V’(— 00 ) = — 00 . Hence by Lemma 6.1 and Lemma 6.2 in the Appendix, we arrive at 
the conclusions. □ 


3 Invariant density 


In this section, for Rq > 1 , we give the existence, uniqueness, global asymptotic 
stability and an explicit formula of the invariant density of the FPE associated with 


SDE (1.3) 


(3.1) 


dp{t, x) 
dt 


~ 0a;]p(t,a;)} + 2 ^ q^2 


^ [x‘^{N — x)‘^p{t,x)). 


Theorem 3.1. If Rq > 1, then there exists a unique invariant probability measure 


for (3.1) which has the density p® with respect to the Lebesgue measure. Moreover, 


the invariant density p^. is globally asymptotically stable in the sense that 


-V 


(3.2) 


lim 

t^OO 


\V{t)q{x) -pI{x)\<1x = 0, V g G Ll{{0,N)), 


where {V{t)}t^Q is the Markov semigroup defined by (3.1) and L\{{0,N)) := {to G 
L^(M) : vj(x)dx = 1, vj(x) = 0 for x ^ N or x ^ 0, and w{x) ^ 0 for x G M}. In 
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addition, the process It has the ergodic properties, i.e., for any u^-integrable function 

K : 

(3.3) ^^0 f / K{QdT = K{y)iol{dy'^ = 1, 

for all Iq G (0, N). Moreover, the unigue invariant density p^^ of the Markov semigroup 
by 


(3.4) 

with 

(3.5) 


pl{x) := CN^ 


^co(Rq-1)-1 
(— 1)+3 


e 


1 _ ~^o{Ro—d 

O — Cq 


(rIY + - 1) r(c„(fl“ -1)), 


^0 — r)-) the gamma function. 


To prove Theorem 3.1, we need to first study the following FPE associated with 


SDE (2.2) 
(3.6) 

du{t,^) ^ 
dt 




2 ^ 


■(/i+ 7 )e^d 


1 + e€ 


u{t,0 


1 

2 ae ' 


Denote by {W(t)}t^o the Markov semigroup (also called stochastic semigroup, see 
Remark 11.8.1 on p.370 in (^) dehned by ( |3T ). 

Now we give the existence, uniqueness and asymptotic stability of an invariant 
density of the Markov semigroup {U{t)}t^Q. 


Theorem 3.2. If Rq > 1, then there exists a unigue invariant probability 
for (3.6) which has the density m* with respect to the Lebesgue measure, 
the invariant density is globally asymptotically stable in the sense 


measure 

Moreover, 


(3.7) 


lim 

t^OO 


\U{t)v{f) 


-<(0|de = 0, VneL^(M), 


where L\{R) := {/ G L^(R) : f(x)dx = 1 and f{x) ^ 0, forx G M}. In addition, 

the process Yt has the ergodic properties, i.e., for any n^^-integrable function F : 


(3.8) 


Pvo 


1 /■* 

lim - / F(Yr)dT 
t^oo t Jo 


F( 7 )<(dp) 


1 , 


for all To ^ Moreover, the unigue invariant density of the Markov semigroup 
{U{t)]t^o is given by 

(3.9) <(0 := 


where C and cq are defined in (3.5|) in Theorem 3.1 


7 





















Proof. Since Yt is conservative and non-degenerate (i.e., -a'^N'^ >0, V ^ ^ 
by (see Chapter III), there exists a unique classical fundamental solution to 


(3.6) (see also p.l53 in (^). Thus by Theorem 6.3 in the Appendix, there exists a 

for all f > 0, f E 


generalized solution u{t,f) E 
density uq E L\ 


provided that the initial 


For simplicity, we denote by throughout this proof. We hrst give the 
existence, uniqueness and global asymptotic stability of Let Id(^) ;= -1- 


with ao = 


(3N — p — 7 1 


-. Notice that ao > 0 for > 1. It is straightforward to 


^2jY2 

verify that C is a Lyapunov function dehned in the Appendix. Let 5 = 2, ai = ci = 1, 

1 " " 

(3N — /i — 7 — 

-. Then there exists sufficiently large a 2 > 0 such 


C2 = 




that both inequalities in (6.5) in Theorem 6.4 hold. Hence by Theorem 6.4 in the 
Appendix, we have the uniqueness and global asymptotic stability of the invariant 
density u^. Next, we verify that given by (3.9) is an invariant density. From (3.6), 
it suffices to show that 
(3.10) 


(^/ 3 N - p - 7 - - (h + 7)e^ + 




1 -I- e? 




! dV(0 

de 


= 0 . 


In fact, is a solution of the ODE 


(3.n) I ^ 


(^f 3 N - /i - 7 - - (h + 7)e^ + 


0" 


■N^e^ 


1 + 




■M, 


which implies u® solves (3.10). Note that 0 < F(co(i?f — 1)) < 00 for > 1, and 
thus C dehned in (3.5) is hnite. Hence is an invariant density. □ 


By Theorem 3.2 and Theorem 6.7 in the Appendix, we prove Theorem 3.1 


4 Persistence Threshold Theorem 


In this section, for Rq > 1, we dehne a persistence basic reproduction number, disease- 
extinction with large probability and disease-persistence with large probability in term 
of the prohle of the invariant density given by (3.4) and establish a persistence 
threshold theorem. 


Let Rq := 


13 N - 

/i + 7 


be the persistence basic reproduction number. Model (1.3) 


is disease-persistent with large probability if the invariant density peaks at a positive 






























number (or equivalently, the mode of the stationary distribution is positive); other¬ 
wise, if the invariant density peaks at zero (or equivalently, the mode of the stationary 


distribution is zero), then model (1.3) is disease-extinct with large probability. 


In the following, we investigate the prohle of the invariant density p*. 

Theorem 4.1 (Profile of Invariant Density). Assume Rq > 1. Then 


(1). Suppose Rq < 1. Then 
Moreover, 


hmp^(x) = 0. 

x^N 


limp*(a;) = oo. 

a::4'0 


(a) if Rq ^ \ then pI is strictly decreasing in (0,A^); 


Co 


(b) if Rq > then p^^ is strictly decreasing in (0,/_), (/+,A^), and 


Co 


increasing in (/_,/+), where 


(4.1) 




(4 - <co) ± ^(4 + <co)2-16co 


and Cq is defined in (3.5) in Theorem 3.1 


(2). Suppose Rq = 1. Then 


limp* (x) = C/N, 

xlO 


where C is defined in (3.5) in Theorem 3.1. Moreover, 


(a) if Cq ^ 4, then p* is strictly decreasing in (0,iV); 

(b) if Cq < 4, then p* is strictly increasing in ^0, ^1 — ^ ) iV ) and decreasing 
in I ( l-^^iV,ivV 


(3). Suppose Rq > 1, then 


limp* (x) = 0. 


Moreover, p* is strictly increasing in (0,/*(cr)) and decreasing in {R{a),N), 
where 


(4.2) 


Ua) = ^ 


(4 - <co) + ^(4 + i?^co)2-16co 


In particular, the invariant density peaks at R{a). 
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Proof. Note that 


pi{x) = 

where h^{x) = Co{Rq — 1) In a; — [co{Rq — 1) + 4] ln(A^ — x) — Cq- 


X 


, . The limit 

N — X 

lima;-|.^rp^(a;) = 0 follows from lima-^jv h„{x) = —oo. Straightforward calculations show 
that 

_ -4^' + (4 - R^Co)Nx + « - l)coiV2 

^a\^) ~ 


{N — x)^x 


Next, we prove this theorem case by case. 
(1). For Rq < 1, we have 


lim hcrix) = oo, 
x4.o 


which implies 
By Rq = Rq + Cq ^ > 1, we have 


limp*(a;) = oo. 


< 1 _ i < <, 

Co Co 

Hence if Rq < 1 and Rq ^ we have 

Co 


(4 ~ Rq Co) + 16(i?g — l)co ^ 0. 


Thus 


h'^{x) ^ 0 in (0, N), 

which indicates that is strictly decreasing in (0,iV). 

If Rq < 1 and Rq > then 

Co 


h'^{x) < 0 in (0,/_), (/+,A^), and h'^{x) > 0 in (/_,/+), 


where I± is defined in ( 4.1[ ). Hence is strictly decreasing in (0,/_), (/+, A^), 
and increasing in (J_, /+). 


(2). For R^ = 1, 


which implies 


lim ho-(x) = —41niV, 

a;-)-0+ 

limp* (x) = C/N. 

X.10 
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If Rq = 1 and Co ^ 4, then 

h'^{x) < 0 in (0, A^), 


which indicates that is strictly decreasing in (0,iV). 
If Rq = 1 and cq < 4, then 


h'^{x) > 0 in 



and h'^{x) < 0 in 




Hence is strictly increasing in ^0, 1 strictly decreasing in 

(3). For Rq > 1, 

h'^{x) > 0 in (0,4 (ct)), and h'^{x) < 0 in {R{a),N), 

where R{a) is dehned in ( |4.2 ). Hence p^ is strictly increasing in (0, R{cr)) while 
decreasing in {R{a),N), and peaks at /*(a). Moreover, \mixiQh(j{x) = —oo 
implies 

= 0 . 

x‘4.0 


Now we complete the proof. 


□ 


Fignre [T] well illnstrates resnlts in Theorem 4T 


Althongh the mean and variance of the stationary distribntion are given in , for 
sake of integrity, we restate them as a corollary of Theorem 4.1 and give an alternative 
proof by considering the FPE (3.1) directly. 


Corollary 4.2 (Mean and Variance) 

variance of the invariant density pi are given by 


Assume that Rq > 1. Then the mean and 


(4.3) 
and 

(4.4) 


E|/l = r(a) 


Var|/i = (r-r(<T))r(<T), 


where I*{a) = 1 — 


Rlf + 1 


Ro 


<N and I* =[ 1- ]N. 

Rn 
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Figure 1: Profiles of invariant density in different cases. 


Proof. It follows from (3.1) that 


(4.5) x[PN — /i — 7 — /3x]pI{x) — (x)) = 0. 

^ LltZ/ 


Note from (3.4) and Theorem 4.1 

(4.6) 


hmp^(x) = lima;j 9 ®(x) = 0 . 

x^N a:4,0 


Integrating on both sides of (4.5) from 0 to iV, we have 

pN 


x[(3N — yu — 7 — / 3 x]j 9 ® (x)dx = 0, 


'0 


which implies 

(4.7) 


-(V 


-(V 


(/SiV —/i — 7 ) / xp^^{x)dx =/3 / x pl.{x)dx. 
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Multiplying — on both sides of (4.5) and integrating from 0 to N, we have 

X — 


-v 


(4.8) j [^N-^-j-/3x]pldx--a^ 


N 




x{N— xYp‘1{x) + / {N—xYp^Y^)dx 


= 0. 


Combining (4.6)-(4.8), we have 


-V 


E[J] = / xp’ldx = r{cr) 

Jo 


and 


/•AT / \ 2 

Var[J] = / x'^p^dx — ^ / xp^dxj = (/* — /*(a))/*(a). 


□ 


From Theorem 4T it follows the persistence threshold theorem with respect to 
disease-extinction with large probability and disease-persistence with large probabil¬ 
ity. 

Theorem 4.3 (Persistence Threshold Theorem). Assume Rq > 1. 


Co 


(i) If Rq < 1 or Rq = 1 ^ , then (1.3) is disease-extinct with large probability. 


Co 


ii) If Rq > 1 or Rq = 1 > —, then (1.3) is disease-persistent with large probability. 


For Rq > 1, even there exists a stationary distribution as proved in the chance 
for the disease to go extinct may still be very large (or the chance for the disease to 


persist may still be tiny) if Rq < 1 or Rq = 1 ^ The disease will persist with large 

probability only when Rq > 1 or Rq = 1 > ^. Such result interprets differently from 

deterministic epidemic models in that for Rff > 1, the disease is endemic and there 
is no chance for the disease to vanish. However, in the stochastic epidemic model, 
even when the stochastic basic reproduction number Rq > 1, there is still chance 


for the disease to go extinct. Model (1.2) is also different from stochastic epidemic 


models in term of a hnite Markov chain, which concludes that disease goes extinct 
with probability one regardless of the basic reproduction number [^. 


From Theorem 4.3, for Rq > 1, the mode of the stationary distribution is located 


at /*(a). We dehne the stochastic disease prevalence for the SDF SIS model (1.2) by 


cr 


N 


4 — Rq Cq d- y^(4 -|- RqCqY — 16Co 
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Figure 2: Disease prevalence. The dashed line and the solid line stand for the deterministic 
disease prevalence curve and the stochastic disease prevalence curve, respectively. The two 
curves intersect at Rq = 2 . 


a number independent of the total population size N. Recall that for the deterministic 


r I f3N 

model, the disease prevalence is — = 1 — — 7 ^, where Rn = - 

N ° /i + 7 


is the deterministic 


basic reproduction number and I* 
two numbers for > 1 . 



In the following, we compare the 


Theorem 4.4 (Disease Prevalence). Assume Rq > 1. 


( 1 ) If Rq < 2 , then 

R{(j) I* 

N N' 

i.e., the stochastic prevalence is smaller than the deterministic prevalence. 

(2) If Rff = 2, then 

R{a) ^ I* 

N TV’ 

i.e., the stochastic prevalence equals the deterministic prevalence. 

(3) If Rq > 2, then 

i.W) /• 

N N’ 

i.e., the stochastic prevalence is bigger than the deterministic prevalence. 


Proof. Equality cqRq = cqRq + 2 implies 
hia) 1 " 


N 8 


6 — CqRq + {cqRq + 2)^ — 16Co 


IJ(t) I* 

Notice that — 7 ^ < — is equivalent to 


N N 


(4.9) 


(cqR^ + 2)2 — 16co < 2 — + coi?| 
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Since Rq = Rq — 2cq ^ > 1, we have 


2 - + cqRq 

Kn 


^■0 

co ( i ??)2 + 2<-8 


D 

0 

D 




> 


Co«)2 + 2<-4co«-1) 


RF 


co«-2)2 + 2i?, 

RF 


D 


> 0 . 


Hence inequality (4.9) is equivalent to 

64 


-16co < 


W )2 Rg 


16 

y^(2 + coO. 


Note that 


we have 


■16cn 


64 16.. „n. 32(i?^-2) 


(Rgy Rg 


+ ^(2 + CqRq ) — 


(RFy 


R{a) I* ^ jjD ^ n 

— < 2 . 


□ 


Figure demonstrates results in Theorem 4A. From Theorem 4H, we see that 
though the randomness in the transmission coefficient stabilizes the disease-free equi¬ 
librium for Rq < 1 and still reduces severity of disease for 1 < Rq < Rq < 2, it 
enhances severity for Rq 


IF > 2 . 


In other words, when the deterministic basic reproduction number is small, stochas- 
ticity in the transmission rate decreases the prevalence. However, when the deter¬ 
ministic basic reproduction number is big, such stochasticity in fact increases the 
prevalence. 

Persistence result in Theorem 14.31 is different from that stated in Theorem 5.1 
in [^, which claims that almost all the sample paths will fluct uate around the level 
J*(cr) in the sense of inequalities (1.4). Whereas Theorem |4.3| states that it is most 
likely that the disease persists at the level /*(ct). In other words, the most possible 
number of infected individuals in the total population at the stationary distribution 
is R{a). 

Now we compare R{a) and /*(ct). 


Theorem 4.5. Assume RK > 1. Then 


15 















in S 2 

(a) 4(a) > 4(a) if yCg ^ < 4^ < - + -Cg ^ 


(b) 4(a) = h{a) z/4^ = ^ + < J; 


(c) h{a) < h{a) if 1 + 2 cq^ < Rj^ ^ yc^ b 


Proof. Recall that Rq = Rq — 2cq Rewrite R{a) and 4(o'), we have 


and 


W = f 


w = f 


— 4co — CqRq + 2 


D 


2 + cqRq \ 3 cq-Rq 

—^- I — Co + - — 


Hence 

4(a) > R{a) 
Note that 

This implies 

4(a) > IJa) 


2 - 4cn > 


3co4? - 2 


3coR^ - 2 3coR^ + 4 


+ 


> 0 . 


2 + co-R^ ^ 




<t»3co < 
■t»3co < 


3co4? - 2 



Co(3co4? - 2)/2 


2 + co-R( 


D 


'2 + coR{fy 

- I —Co 




'2 + Co4?V 3cn4?-10 


- 


Co < 


4 

-o-fi-o 
12 


10 _1 3 2 

2 ^0 < -Ro < 2 + gCo . 

Combining 4^ = 4^ — 2co ^ > 1, we arrive at the conclusions. 


□ 
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Figure 3: Disease prevalence. In (a), the surface in mixed color stands for I^{a)/N 
and the surface in black stands for I^{a)/N. In (b), (c) and (d), the dashed line 
and the solid line stand for I^{a)/N and h{a)/N, respectively. Two curves almost 
coincide when cq is sufficiently large. 


Results in Theorem |4.5| are illustrated in Figure We see that when cq are 
sufficiently large (or equivalently, a is sufficiently small, ), and almost 

coincide. In fact, it is easy to verify that 


lim 4(cr) = lim4((T) = I*, 

(J— )-0 cr^O 


which is a corollary of Theorem 5.1 in the next section. 


5 Limit Stochastic Threshold Dynamics 


In this section, we further investigate the asymptotic dynamics of (1.3) and try to 
establish a sharp connection between (1.3) and its deterministic counterpart (1.1) in 
term of the limit of the invariant density as a —)■ 0. 


Theorem 5.1 (Limit Stochastic Threshold Theorem). 


(a) . //< ^ 1, then the number of susceptibles I at DEF converges to 0 in probability 

as cr —0. 

(b) . If Rq > 1, then the number of susceptibles I atp^. converges to I* in probability 

as cr —)■ 0. 
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Proof. Note for Rq ^ 1, 


Rq < 1, V cr > 0. 


Conclusion (a) trivially follows from Theorem 2.2 for the random variable / at DFE 
is identically zero. 

In the following, we only prove (b). 

It suffices to show that for small £ > 0, 

pl* -\-£ 

/ p^^{x)dx —)■ 1, as (T —)■ 0. 


(5.1) 

Since 

(5.2) 

for sufficiently small a, 

(5.3) 

This implies 


' I*-e 


lim R{a) = /*, 

( 7 - 1-0 




(5.4) I* — e < R{a) —I* + e > R{a) + 


Since Rq > 1 and Rq = Rq -—, we know that Rq > 1 for all sufficiently small 


(T^. By Theorem 4.1, p* is increasing in (0, and decreasing in {R{a), N). Thus 

/•/*(cr)-£/2 

/ pl{x)dx 

Jo 

<{R{a) - e/2)pl{h{a) - e/2) 

= (4(a) - e/2)pl{R{a) - £/4)e^-C*C)-^/2)-/i.T*C)-^/4) 


R{a) ^/‘^ ^ha(U(a)-el2)-h„(h(a)-el‘i) 

pl{x)dx. 


< ' ' —I p;(a;)dx 

J/*((7)-£/4 

R{a) - £/‘J h^(I4(7)-el2)-h^(I,{a)- ■ ' 

e/4 


-£/4) 


'/*(o-)-e/2 


Similarly, 


[ pl{x)dx < - -- ^J^^h^{I4P+e/2)-h^{I,{a)+e/4) ^ ^ pl{x)dx. 

'h{a)+e/2 ^/4 Jh(a)-e/2 
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Hence by (5.4), we see 


'■h{o)-e/2 


rN 


r‘/t(cr)+£/2 


1 = 


p’^{x)dx + / p%{x)dx + 


^ 1=^ [a )+£/2 


fh{a)-e/2 


p® (a;)dx 


< 


^ 4(0-) ^hM,M-el2)-h^(h(a)-6l4) 

e/4 


< 


^ N U{a) ^h^(I,(a)+e/2)-h^(h{o)+elA) 
e/4 

^ _j_ ^/‘^ ^ha{h{a)-e/2)-haiIt{a)-£/A) 

e/4 


^ 7 * (cr )+£:/2 


^7*((t)— e/2 


p® (a;)da; 


+ 


^ I*{o') ^/‘^ ha-ih(a)+e/2)-h^ih(a-)+e/4:) 

e/4 


-/*+£ 


p® (a;)da;. 




In the following, we only need to show that 


£/ 4 ) 


e/4 


I^^c) ^/‘^ ha-(h((T)+£/2)—ha-(h{a-)+e/4:) 

e/4 


—y 0, as (j —y 0. 


By (5.2), it suffices to prove that 

(5.5) h„{h{a) — e/2) — h„{h{a) — e/4) —)■ —oo 

and 


(5.6) h^{I^,{a) + e/2) - ha{h{a) + e/4) -)■ -cx). 

Straightforward calculations give 

Since 

(5.7) Co —)■ oo, as cr —)■ 0, 
which implies 

hnih'/(4(cr)) = -oo. 

(T—^0 

Hence for small e > 0 and sufficiently small a > 0, h'^ is decreasing in (/*(cr) 
e/2,4((T) +e/2). Thus 

K{h{a) - e/2) - K{h{a) - e/4) < -|h/(4(cr) - e/4). 
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Recall that 


Kix) = 


-4x^ + (4 - R^co)Nx + « - l)coiV^ 
(iV — xYx 


Substituting = 0, by (5.7), we have 


KiW 


£ X —+ Sslif ((t) + eN(^Rq cq — 4) 
4^ ^ 4[N - (4(a) - £/4)]2(R(a) - e/4) 


—)■ oo, as (T —)• 0. 


Now (5.5) is proved. Using similar arguments, we can prove (5.6). 


□ 


Recall that for deterministic SIS model (1.1), we have the threshold theorem: the 
disease-free equilibrium Pq = (iV, 0) is globally asymptotically stable if Rq ^ 1, and 
a unique endemic equilibrium P* = {N — I*, I*) is globally asymptotically stable if 


Rq > 1. Theorem 5.1 establishes a sharp link between the SDE SIS model and ODE 


SIS model in term of their threshold dynamics as the randomness in the transmission 
coefficient vanishes. 


6 Summary 


In this paper, we further study the global dynamics of an SDE SIS epidemic model 
proposed in |^. Using Feller’s test for explosions of solutions to one dimensional SDE, 
we establish a stochastic threshold theorem and thus prove the conjecture proposed in 
[^. By studying the FPE associated with the SDE, we prove the existence, uniqueness 
and global asymptotic stability of the invariant density of the FPE. Using the explicit 
formula for the invariant density, we dehne the persistence basic reproduction number 
and give the persistence threshold theorem in term of the invariant density. By 
comparing the stochastic disease prevalence with the deterministic disease prevalence, 
we discover that the stochastic prevalence is bigger than the deterministic prevalence 
if the deterministic basic reproduction number Rq > 2. This shows that the disease 
is more severe than predicted by the deterministic model if there is randomness in 
the transmission coefficient when the basic reproduction number is large enough. 
Finally we investigate the asymptotic dynamics of the SDE SIS model and establish a 
connection with the dynamics of the deterministic SIS model as the noise vanishes. We 
expect such global stochastic threshold theorem, persistence threshold theorem and 
limit stochastic threshold theorem discovered in this simple SIS type model to exist 
in more complicated epidemic models, for instance the SDE SIR model [^[^. We 
also expect that the approach used in this paper applies to other (high dimensional) 
SDE biological models, by considering the associated FPE and its invariant density, 
even though the dynamics may be more complex and there may not exist an explicit 


formula like (3.4) for the invariant density. Nevertheless, Feller’s test may not work 


for Stochastic Threshold Theorem (like Theorem 2.2) if the model can not be reduced 
to a one-dimensional SDE. 
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Appendix 


First, we present the monotonicity of function g dehned in (2.1[) in Section 2. 


Lemma 6.1. Thefunctiong definedin (2.1) is strictly increasing in {0, N). Moreover, 


limgf(f) = —oo, and limgf(7) = oo. 

ao 


Hence, its inverse function g~^ exists and is strictly increasing in (—oo, oo). In addi¬ 
tion, 

lim g~^{y) = 0, and lim g~^{y) = N. 

y^—oo y—^oo 
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Now we give some preliminary results on Feller’s test of explosions and FPEs 


10,15 


The following is a result on Feller’s test of explosions [7, 11 . For a complete 
statement of Feller’s test of explosions, we refer the reader to Proposition 5.22 on p. 


345 in 11 


Consider the following one-dimensional SDE 


( 6 . 1 ) 


dX(f) =b{X{t))dt + a{X{t))dB{t), 

X(0) =Xo 


for some Xq G M. Assume that 
(HI) {a{x)f >0, Vx e M. 

(H2) VxgM, 3£>0 such that / ^ d,^ < oo. 


' x—e [a, 


■«))■ 


Dehne the scale function by 


ip{x) = / exp 

Jo 


2b{r) 


rdr d^, a: G M. 


Jo (a(r))' 

Now we give Feller’s test [II] . 

Lemma 6.2. Assume that (HI) and (H2) hold. Let Xq G M. 

(1) //'0(—oo) > —oo and 'ip{oo) = oo, then 

p| lim X{t) = — oo| = 1. 

L t^oo J 

(2) //'0(—oo) = —oo and 'ip{oo) = oo, then 

p| sup X{t) = oo| = p| inf X{t) = — oo 

In particular, the process Xt is recurrent: for every G M, 

P{X(f) = e : 3 f G [0,oo)} = 1. 


= 1 . 


Next, we present some basic results on FPE. Consider the following initial value 
Cauchy problem: 

f)v f) 1 / 9 \ 

(6.2) .- = -^^b{x)v) + -—(^{a{x)) f>0, xgM, 
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(6.3) 


f(0,a;) = Vo{x), a: G M. 


Eq.(6.2) is the Fokker-Planck associated with (6.1). It is well-known that there is a 


unique generalized solution to (6.2) provided that there exists a unique classical funda¬ 


mental solution to (6.2) ( for the dehnitions of fundamental solution and generalized 
solution, see, for instance, p.365 and p.368 in p^). 


In the following, we assume there exists a unique classical fundamental solution 


to (6.2) and thus 


Theorem 6.3. 


v{t, x) to (6.2). 


15 


For every vq G there exists a unique generalized solution 


Eq. (6.2) dehnes the Markov semigroup (also called stochastic semigroup in [15 


{V(t)}t^o of operators on EC,_(M) by 

(6.4) V{t)vo{x) = v{t,x), 

and the set of all densities E)^(]R) is invariant under V{t) for each f > 0. 


For more preliminaries on FPEs, we refer the reader to 15 


Now we state a known sufficient condition given on p.742 in [1^ for the existence, 
uniqueness and global asymptotic stability of an invariant density (see also Theorem 2 
and Theorem 4 on p.l54 in and Theorem 11.9.1 on p.372 in [Is]). 

We dehne a Lyapunov function V : M —)■ M as a function with the following 


properties 17 : 


(PI) V{x) > 0,V a; G M; 

(P2) 1^13,1^00^(0;) = cx); 

(P3) V{x), |E'(a;)| ^ 5ie^'^'^^\'i x G M, for some positive constants ^i, 52- 


Theorem 6.4. Suppose there exists a Lyapunov function V and positive constants 
Cl, C 2 , 5, cxi and a 2 such that 


— cti -t- Cilxl”^ ^ Vix), 
b{x)V'{x) -F ^a^{x)V"{x) ^ -C 2 V{x) -f 02 - 


Then there exists a unique invariant probability measure p for (6.1) which has the 
density fo with respect to the Lebesgue measure. Moreover, the invariant density /o 
is globally asymptotically stable in the sense that 


( 6 . 6 ) 


lim 

t^OO 



|V(t)/(a;) - /o(a;)|dx 


0, V / G L\{R). 
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In addition, the process Xt has the ergodic properties, i.e., for any p-integrable func¬ 
tion F : 


(6.7) P.( i ^ F(X,)dr = J F{y)p{dy)) = 1, 

for all Xo = X G M. 


Remark 6.5. In 15 , regularity oi coefficients are required in Theorem 11.9.1. However, 
as pointed out on p. 365 in 15 , this regularity is simply for the existence and 
uniqueness of solutions to (6.2) and (6.3). In other words, Theorem 11.9.1 on p.372 
in 15 without regularity condition still holds if the existence and uniqueness of 


solutions to (6.2) and (6.3) are provided. 


In the following, we state the relationship between the solution u{t,^) to (3.1) and 
the solution p(t, x) to (3.6). 

Theorem 6.6. 


p{t,x) = g\x)u{t, g{x)) = 


N 


u 


x{N — x) 

provided that the initial densities satisfy 


a; ) ’ allt>0, x G (0, N) 


Po{x) = , ^oflog 


x{N — x) 

In particular, the invariant densities satisfy 

. N 


N — X 


PAX) = 


ul log 


X 


x{N — x) —X/ 

Proof Let Fi(t,x) = J^p{t,y)dy and F2{t,f) = J^^u{t,p)dp. Then 
Fi{t,x) = F2{t,g{x)), V t > 0, X G (0,X). 
Hence for alH > 0, x G (0, N), 

= -^{t,9{x))g (x). 


i.e.. 


Pa{t,x) = g'{x)u{t,g{x)). 


□ 


Theorem 6.7. 
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(1) K is an invariant measure for (2.2) if and only if v = k o g is an invariant 
measure for (1.3). 


(2) u^. is asymptotieally stable if and only if p^. is asymptotically stable. 

(3) Yt is ergodic if and only if It is ergodic. 


Proof. We only prove case (2). Suppose is asymptotically stable, i.e., 


lim 

)-oo 


|W(t)/o(0-<(0|de = 0, V/oGL^ 


V Wo e L^(0,A^)), let /o(0 = [g HO)) then it is easy to verify that 

/o € Moreover, we can show that 



\V{t)wo{x) 


p"^{x)\dx 


/ OO 

|w(t)/o(0-<(0|d0 

-OO 


which implies that 


rN 


lim / \V{t)wo{x) — pl.{x)\dx = 0, 

t^oo Jq 

i.e., pI- is asymptotically stable. Conversely, we can also prove that is asymptoti¬ 
cally stable provided that is asymptotically stable. □ 


26 




